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The d-wave vortex lattice state is studied within the 
framework of Bogoliubov-de Gennes (BdG) mean field the- 
ory. We allow antiferromagnetic (AFM) order to develop self- 
consistently along with d-wave singlet superconducting (dSC) 
order in response to an external magnetic field that generates 
vortices. The resulting AFM order has strong peaks at the 
vortex centers, and changes sign, creating domain walls along 
lines where V x j s w 0. The length scale for decay of this AFM 
order is found to be much larger than the bare d-wave coher- 
ence length, £. Coexistence of dSC and AFM order in this 
system is shown to induce 7r-triplet superconducting order. 
Competition between different orders is found to suppress the 
local density of states at the vortex center and comparison to 
recent experimental findings is discussed. 



I. INTRODUCTION 

Many of the anomalous properties of high T c supercon- 
ductors (SC) are believed to be due to the competition 
between different kinds of orders. The electronic struc- 
ture of vortices in the cuprates is one striking example 
of behavior which is not understood and which may be 
influenced by competing or coexisting order parameters. 
The electronic structure of s-wave vortices has been stud- 
ied extensively by both theoretical |l],^] and experimental 
H methods which established the existence of localized 
vortex states with very low energy. The situation is more 
complex for the case of d-wave SC due to the presence 
of nodes in the energy gap. Earlier theoretical results of 
Wang and MacDonald performed within the frame- 
work of Bogliubov-de Gennes (BdG) calculations, pre- 
dicted a strong and broad zero bias conductance peak 
(ZBCP) in the local density of states (LDOS) at the vor- 
tex core. This ZBCP was shown || to display a four-fold 
symmetric "star-shape" , oriented towards the gap nodes 
and decaying as a power law with distance from the vor- 
tex core. In contrast, experiments present a completely 
different picture. Tunneling studies of both YBCO ||] 
and BSCCO § show no evidence of a ZBCP. Further- 
more, structure is observed in the low energy differential 
conductance data which has been interpreted as the ex- 
istence of low energy quasiparticle states at ±5.5 meV 
(for YBCO) and at ±7 meV (for BSCCO). These states 
persist to long distances - a size much larger than that 
of the vortex core. 

There have been a number of different proposals for 



reconciling this contradiction between theoretical pre- 
dictions and experimental findings, including effects of 
strong correlations and competing order parameters. The 
role of strong correlations was emphasized by Franz and 
Tesanovic M who argued that the correct description of 
the normal state of high T c cuprates is a doped Mott 
insulator and not a simple metal (as in BCS theory). 
Hence, at the vortex core, where superconductivity is 
suppressed, a gapped insulator should be expected. In 
this case, the pseudogap behavior of the tunneling re- 
sults at the core center is a result of spin-charge separa- 
tion. To address the correlation effects on the structure 
of an isolated vortex, Han and Lee Jl(j studied the t — J 
model using resonant valence bond mean field theory. 
Their main result is the prediction of two types of vortex 
core - insulating or metallic, depending on the model pa- 
rameters. They also found evidence of staggered current 
patterns Q in the vortex core. The pairing order param- 
eter near the vortex core was found to have d x 2_ y 2 +id xy 
symmetry [[l0 11 . A different argument for the absence 



of the ZBCP invokes the structure of the c-axis tunnel- 
ing matrix element Jl^ ]. In a recent study, Berthod and 
Giovannini |]l3| showed that with increasing correlation 
strength the BCS core states become suppressed and the 
spectra inside and outside the vortex core become simi- 
lar. 

An alternative scenario has been presented in which 
antiferromagnetic (AFM) order is generated near the vor- 
tex core. Arovas and collaborators Jl4|] first studied this 
problem within the framework of SO(5) theory and pre- 
dicted that the vortex core becomes antiferromagnetic, 
and hence insulating, at low temperature and relatively 
low doping. Subsequently, a BdG calculation of SO (5) 
vortices ]15| found that the AFM order, although itself 
gapless within this model, suppressed the ZBCP. Recent 
works of Sachdev and collaborators |l6| argue that the 
system undergoes a quantum phase transition from a 
pure d-wave superconducting (dSC) phase to a mixture 
of dSC and spin density wave (SDW) phase as a function 
of the external magnetic field. In this theory the dynamic 
spin fluctuations in dSC are pinned by the vortex cores. 
More recently Zhu and Ting [l7| have studied a d-wave 
vortex lattice, allowing for AFM order, using BdG mean 
field theory. They find AFM order which alternates in 
sign from one vortex to the next, and a splitting of the 
ZBCP due to the AFM order. 

There have been a large number of recent experiments 
which suggest the emergence of AFM or SDW order on 
top of the dSC order in the cuprates in an external mag- 
netic field. Using inelastic neutron scattering, Lake et al. 
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|L8| found that the magnetic field induces a large increase 
in the spectral density of low energy spin fluctuations in 
LSCO near optimal doping. This demonstrates the ex- 
istence of a fluctuating magnetisation density with a pe- 
riodicity of 8 lattice spacings along the direction of the 
Cu-0 bonds. Scanning tunneling microscopy (STM) on 
slightly overdoped BSCCO in the mixed state shows 
a charge modulation in the LDOS data that has a peri- 
odicity of 4 lattice spacings - a result consistent with the 
neutron scattering study since the periodicity of charge 
modulation is one half of that of the SDW modulation. 
Field induced enhancement of the Bragg peak intensity 
has also been observed by Khaykovich et al. [Q, in an 
elastic neutron scattering study of La2Cu04+ y , indicat- 
ing a field induced static AFM order. High-field nuclear 
magnetic resonance (NMR) studies of YBCO [^fj have 
found evidence for strong antiferromagnetic fluctuations 
inside the core. Muon spin rotation measurements ]22] ] 
have also been interpreted as evidence of static mag- 
netism in the vortex core of Ortho-11 YBa2Cu3 06.5- Un- 
derstanding the charge and spin modulation in the vor- 
tex lattice state is currently under intensive theoretical 
investigation fH-§j|. 

Motivated by the discrepancy between theoretical pre- 
dictions and experimental observations, we study the 
structure and spectrum of a d-wave vortex lattice us- 
ing a BdG mean field theory which allows both static 
AFM order and dSC order in a self-consistent manner. 
Our main results are: (1) In the lowest energy configura- 
tion, AFM order has the same periodicity as the vortex 
lattice, in contrast to earlier studies |I7|, although there 
are antiphase domain walls in the AFM order within the 
unit cell; (2) triplet SC order is induced in the presence 
of coexisting AFM and dSC, which approximately sat- 
isfies the constraint predicted by Zhang |2(|; (3) due to 
the presence of AFM order, the LDOS exhibits structure 
which persists well beyond the vortex core size. 

The organization of the paper is as follows. In Sec. II, 
we describe the class of models we study, and we also 
elaborate on the BdG method. Sec. Ill discusses our 
results on the spatial structures of different order param- 
eters in the mixed state, and Sec. IV discusses the LDOS 
results. Finally, we present conclusions in Sec. IV. 



II. MODEL AND METHOD 

We describe a 2D d-wave SC in an external magnetic 
field by the following Hamiltonian : 



Tt = JC + Hint + Wp, 

where the kinetic energy, 

K = ~ t E (e^4 a c ja +h.c.), 



(1) 



(2) 



<ij>,a 



describes electrons, with spin a at site i created by c\ a , 
hopping between nearest-neighbors (ij) on a square lat- 



tice. The Peierl's phase factor 0y = (7r/$ ) J ' * A(r) -dr, 
with the superconducting flux quantum $o = hc/2e, de- 
notes the presence of an external magnetic field produced 
by the vector potential A(r). The interaction term in 
Eq. ©, 



J 



<ij> 



Si ■ s , 
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(3) 



is chosen to generate both antiferromagnetic and d-wave 
superconducting orders at the mean field level. Here the 
spin operator is 5f = c.Lc^c^, where the a a are Pauli 

matrices, and the density operator is n, a = c\ a Ci a with 
n>i = Uii + n il- An ons it e repulsion term of the form 
UJ2i n i1 n il can be added to Hint to ensure the absence 
of any onsite s-wave pairing, however, we found that it 
is not necessary because such induced s-wave pairing is 
vanishingly small for our parameter regime. The chemi- 
cal potential part of the Hamiltonian, TL^ = — fi^j rii, is 
used to fix the total density of the system to the desired 
value. 

It is commonly believed that this simple model, in the 
limit of no double occupancy of any site, is adequate 
to describe the strongly correlated cuprates at low tem- 
peratures. A mean field treatment of this model, while 
missing the strong correlation effects, may still usefully 
describe the interplay between static antiferromagnetic 
order and d-wave superconductivity ]27|] . 

A mean field decomposition of the above model leads 
to the following effective Hamiltonian 



H 



eff 



E(*+^ 



y~]\fH- Vi + ae lQri m i+ s J 

i.a \ 8 / 

E( A '[ c U c a + c lT c U] + h.c 



(4) 



where rf a is the Fock shift that renormalizes the hopping 
amplitude t and rji is the Hartree shift that renormalizes 
the chemical potential fj,. Here S = ±x, ±y, the four 
nearest neighbours of any site, and <j>f is same as 0y- 
with j = i + S. In the above equation the local dSC (A*) 
and AFM (rrii) order parameters, and the Hartree (rji) 
and Fock {r io ) shifts satisfy the following self-consistency 
conditions: 



mi = - ((n iT ) - (n a ))e lQ - r * 

T ia = 7r( c i+,5,-o- C i,-o~) 



(5) 
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where Q = (71", 7r) is the antiferromagnetic wave vector. 
The Hamiltonian in Eq. (Q) has been studied for a range 
of model parameters. In the following we will report re- 
sults primarily from our calculations at zero temperature 
(T = 0) with parameters J/t = 1.15 and (n) = 0.875. In 
the absence of an external magnetic field, these parame- 
ters result in a uniform d-wave gap, A max ~ 0.36i (which 
is 8 A in our notation) , with superconducting coherence 
length £0 ~ 4a (where a is the lattice spacing), a value 
consistent with high T c cuprates. 

We have also studied an extended Hubbard model 
where the dSC is generated by a nearest neigh- 
bour attractive interaction of the form: Hint = 
—J'Yl<ij> a a' n iat n ja'- Self consistent calculation on 
such a model will generate a triplet superconducting or- 
der in addition to a singlet dSC order when decomposed 
in a mean field manner. Also, such an extended Hubbard 
interaction does not by itself generate AFM order, which 
has to be incorporated into the model using an addi- 
tional interaction parameter [jl7|. By contrast, both dSC 
and AFM orders are generated from the single J-term in 
Eq. (j^) in our calculations, thereby putting both orders 
on an equal footing. Another problem of the extended 
Hubbard model is that, the strength of the onsite repul- 
sion, U, required to generate AFM order has to be fine- 
tuned to study the interplay of the two competing orders. 
Non-zero AFM order develops only for U > U c , a result 
consistent with earlier findings jl7) . On the other hand a 
strong U would generate strong antiferromagnetism that 
would in turn kill d-wave superconductivity. Except for 
these limitations, if the triplet order is ignored, the mean 
field decomposition of the extended Hubbard model gen- 
erates an effective Hamiltonian similar to that given by 
Eq. ([|), but with different values for the Hartree and Fock 
shifts. We have checked that the qualitative physics de- 
scribed by the extended Hubbard model at the mean field 
level remains similar to that for Eq. (|J), reported in this 
work. 

TL c s can be diagonalized using Bogoliubov tranforma- 
tions, 



{ u «,a(«)7ri,a - a< >a (i)7n,-a} 



leading to the usual BdG equations 



£-aM A 
A* -|* - aM 



^n. — a. 



= E n 



vn.—a 



(6) 



(7) 



where £tin >a (i) = -J2s t i u n,<x( i + S ) ~ fii u n,a{i), 
Mu ni a(i) = (J28 m i+s) etQ ' riu n,a(i) and Au„ iQ (i) = 
X}<s(Af + A^ s )u n ^ a (i + 8). One can numerically diago- 
nalize Eq. (|7]) directly. However, exploiting the magnetic 
translation symmetry of the effective Hamiltonian allows 
one to treat much larger system sizes. This procedure was 
discussed by Wang and MacDonald [Q. However, since 
our model and resulting BdG equations differ somewhat, 
we describe below the block diagonalization procedure 
we use. 



We assume that our system is made up of a periodic 
array of identical magnetic unit cells of size N x x N y . 
For our case we take N x — N y /2, a rectangular mag- 
netic unit cell enclosing two superconducting flux quanta, 
and work in the Landau gauge, A(r) = (0, Hx). In this 
gauge the magnetic translation operators are (r|TR|r') = 
5 r}T +fie~ , ' bR ' : ( lv+Ry \ where the parameter b is defined as 
b = H/$ = 2ir/N x N y . Since [H cS ,T R } = 0, the eigen- 
states of Tn can be used to block diagonalize H e s- The 
transformations that block diagonalize the BdG equa- 
tions in the magnetic wave vector, k, are 



< >a (R) = e lk - R e -' ib (^+A B )A* u nj k ) 



(8) 



where k x — 



with m x = 0, 1, 2, (p — 1) and 



vj, - -jiv Jt ) witn m y = 0)1)2, —,{q - 1)- The eigen- 
value problem for a system of size pN x x qN y reduces to 
p x q eigenvalue matrix equations, each for a system of 
size N x x N y . The block diagonal BdG matrix of Eq. (^) 
becomes, for each k, 



e,(k)-«M(k) A(k) \(ul a (k) 

A*(k) -&(k) - oM(k) J U\-a( k ) 

< Q (k) 



«(k) 



(9) 



Here the operators (in the bulk) are defined as 

fa(kK (k) = -X;{*+^^ +a ,_ tt C jl -„>}. 



Mi uV-, 



,(k) 



Af(k)< Q (k) = 



/ 4 m q.n' 
q,n' 



A(k)<_ Q (k) = 2J2 \Y1 \ w "+<5,-«(k) ■ (10) 

In this 'repeated zone scheme' the self-consistency con- 
ditions, Eq. (|^), can be rewritten as 

Knd) = ^ [K;(q)«?t(q) + <i(qK*(<i)} i 1 - f^i) 

+ |< T (q)^t(q) + ^ T (q)<l(q)} / q „ T ] (11) 

m qn (i) - [K,M)?Ua l<i(q)l 2 /qm 

+ |< T (q)| 2 (l - / qrU ) - |< ; (q)| 2 (l - / q „ T )] (12) 

where j = i + 5 and J\f = p x q, the total number of mag- 
netic unit cells. Note that, for our choice of gauge, the 
matrices £ Q (k),M(k) and A(k) are k-independcnt for the 
bulk of the system, as can be seen from Eq. ([To|). How- 
ever k-dependence enters through the boundary terms 
(See Appendix). 

Starting with an initial guess for all local variables (dSC 
order, AFM order and Hartree and Fock shifts) and for 
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the value of /z, we solve for all eigenvalues and eigenvec- 
tors of the BdG matrix for each k, in terms of which all 
local orders are recalculated using self-consistency condi- 
tions Eq. (fl2|). This process is repeated until a conver- 
gence is obtained in these variables and also in fi which 
fixes the total density to the desired value of 0.875. 

III. RESULTS 

A. Spatial structure of local order parameters 

In this section we discuss the structure of the dSC, 
AFM and charge density orders, as well as of the 7r- triplet 
SC order which is induced in the presence of coexisting 
AFM and dSC orders. 

d-wave Pairing Amplitude: For our lattice BdG cal- 
culation we defined a d-wave pairing amplitude A 5 (i, R) 
following Eq. (E|) that lives on bonds. We can also de- 
fine a gauge invariant site dependent singlet d-wave order 
parameter in our model as, 

Ai = J [A*(r) + A- £ (r) - e - lba: A s (r) - e^A^r)] 

(13) 

where r = (i, R). In the left panel of Fig. (Q) we plot 
the magnitude of Aj on half of a magnetic unit cell. The 
spatial dependence of A, is very similar to its behavior in 
the absence of AFM order (not shown here). The size of 
the vortex core, £, is similar for the two cases, although 
the initial slope inside the core is somewhat smaller in 
the case with AFM order. This is what one would expect 
within Ginzburg-Landau theory |3(^,|3l| . 

The presence of the vortex lattice induces an extended 
s-wave order (A|) around vortex in a self-consistent man- 
ner; defined as: 

Af = [A* (r) + A" £ (r) + e^A^r) + e ibx A-y(r)] /4. 

However, the strength of the induced A s order is vanish- 
ingly small compared to the primary d-wave order, and 
hence it will be ignored in further discussions. 

Antiferromagnetic Order: In the absence of an exter- 
nal magnetic field, the antiferromagnetic order is stabi- 
lized at half filling ((n) — 1) and away from half filling, 
AFM order is quickly suppressed. 
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FIG. 1. Upper panel: Singlet d-wave pairing amplitude Ai 
on a magnetic unit cell of size N = 26 x 26, with J — 1.15t 
and (n) = 0.875. Ai is uniform (and equal to the BCS value) 
away from the vortex core, and falls to zero at the vortex cen- 
ter with the length scale £. 

Lower panel: The grey-scale plot shows the length scale £. 
The dark (light) regions indicate smaller (larger) values for 
Ai. Note that only half of the unit cell, containing one flux 
quantum, is shown here (and in all other figures). The full 
unit cell configuration can be obtained by periodically repeat- 
ing the results along y-direction. 

For our choice of parameters, the minimum energy con- 
figuration is one in which AFM order is identically zero 
for the uniform system at zero external magnetic field 
(H = 0). However, at low temperatures, in the presence 
of a magnetic field, AFM order develops self-consistently 
in the vicinity of the vortex cores where superconductiv- 
ity is destroyed. The structure of the AFM order is very 
different from that predicted earlier |17| ]. We have con- 
firmed that the structure of our order parameter, which 
changes sign within the unit cell as is discussed below, 
minimizes the ground state energy, E g , the expectation 
value of the Hamiltonian of Eq. ([j]), evaluated in the 
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mean field ground state. 
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FIG. 2. Upper panel: Staggered magnetization rm on a 
magnetic unit cell of size N — 26 x 26 system, with J = 1.15t 
and (n) = 0.875. rm is strongest in the vortex cores, and 
decays very slowly away from it. 

Lower panel: The grey-scale plot shows the much larger 
length scale for the decay of rm compared to £ (See also 
Fig. |l|). The dark (light) regions indicate larger (smaller) 
values for m;. 

In this regard, we emphasize the role of including the 
Hartree and Fock shift terms. In a variational calcula- 
tion, such as the BdG method, it is possible to achieve 
self-consistency for various different order parameter con- 
figurations corresponding to different -E s 's. The correct 
configuration is the one that minimizes E g . Ignoring the 
Hartree and Fock shifts will give a higher value for E g . 

The order parameter configuration in Ref. p7| , when 
periodically repeated for an array of unit cells along x 
and y directions, produces an AFM order that changes 
sign on alternate unit cells along the ^-direction, and 
has the same sign along cc-direction. Such a "stripe-like" 
AFM order spontaneously breaks the symmetry of the 
x and y directions. Working with a unit cell contain- 
ing 4 vortices, we have verified that such a "stripe-like" 



AFM order does not converge to a self-consistent solu- 
tion within our Hartree-Fock, BdG formalism. On the 
other hand, a "checkerboard" pattern of AFM order, in 
which the sign of AFM order alternates in adjacent unit 
cells along both the x and y directions, leads to a self- 
consistent solution that has slightly higher E g than for 
the periodic configuration reported here. 

Comparison of Fig. (^) to the result for the d-wave 
order parameter shown in Fig. (Q) above, clearly shows 
that the length scale of decay of the AFM order away 
from vortex core is much larger than the vortex core size 
£, as one would expect in the regime where uniform dSC 
is stable in the absence of an applied magnetic field. 




FIG. 3. Contour plot of staggered magnetization m< (left 
panel) on a magnetic unit cell of size N = 26 x 26 system, 
which shows that rrii changes sign creating a domain wall 
along the line on which V x j 3 « 0. The contours correspond 
to m = 0.075, 0.03, 0.004, 0, -0.01 from the center to radially 
outwards. Here, j s is the supercurrent associated with dSC 
order, which is plotted in the right panel. Notice that along 
the line where rm ~ 0, j s arrows do not bend, showing that 
V X j s « along those lines. 

We also notice in Fig. (^|) the surprising result that 
the staggered magnetization changes sign within the vor- 
tex unit cell. Specifically it changes sign along the line 
where the supercurrent j s associated with the dSC or- 
der parameter has a zero curl, i.e. the line along which 
j s changes the direction of its winding. This result is 
true for a wide range of parameters studied (i.e. J, (n) 
and system size or external field strength); as well as for 
the study with the extended Hubbard model described 
in Sec. ||). This puzzling effect can be thought of as the 
result of an effective interaction of the following form: 
p8[ (V x j)/(V xj)jmjmj, where I and J label adjacent 
plaquettes and to/ is a course-grained AFM order associ- 
ated with the plaquette. Alternatively, this implies that 
the gradient term for the AFM order in the G-L free en- 
ergy takes the form |(V — a)m\ 2 , where a is proportional 
to the supercurrent. Note, that this interaction implies 
no preferred sign for the AFM order, but only that it 
changes sign along the line where the supercurrent vor- 
ticity changes sign. We have confirmed that this is also 
the case for our numerical BdG solution, i.e. there is no 
preferred direction for the AFM order. 



5 



Charge Density Order: In the absence of AFM or- 
der the charge density is relatively structureless; (rij) is 
uniform everywhere except very near a vortex. Around 
the vortex it has a very weak dip which changes to a 
weak peak at the vortex center. We believe that such 
weak density fluctuations would not survive long range 
Coulomb screening effects and hence will not have any 
important effects in a more realistic model. However, 
when the AFM order is allowed to develop self consis- 
tently, the charge density reorganizes itself significantly 
to accommodate the AFM order, as shown in Fig. 





FIG. 4. Upper panel: Charge density (m) on a mag- 
netic unit cell of size N = 26 x 26 system, with J = 1.15t and 
(n) = 0.875. (rii) is more or less uniform away from the vortex 
core, but achieves a value very close to half-filling ({rn) — 1) 
at the vortex center. This local half-filling at the vortex core 
favors the AFM order in the core region where dSC order is 
depleted due to external magnetic field. 

Lower panel: The grey-scale plot shows that the distance from 
the vortex center where (n,) differs significantly from its uni- 
form value is slightly larger than f. (See also Fig. |l|) The 
dark (light) regions indicate larger (smaller) values for (m). 



The local electron density remains more or less uni- 
form away from the vortex center, but near the vortex 
core the magnitude of (n^) increases continuously result- 
ing in a density close to half filling at the vortex center. 
For the t — J model at T — 0, AFM order becomes the 
dominant order at half filling in the absence of the mag- 
netic field. The charge density thus organizes itself to 
favor the AFM order strongly at the vortex core by ap- 
proaching unity. Again, our result is different from that 
in Ref. p7| , where (rij) differs from the uniform bulk value 
only within one lattice spacing of the vortex center. Full 
self-consistency, including the Hartree and Fock terms, is 
required to reliably study charge inhomogeneities. How- 
ever, one does need to keep in mind that the long-range 
Coloumb interaction in a three-dimensional material can 
further alter the charge distribution in ways which are 
beyond the scope of the present study. 

Other self-consistent local variables such as the Hartree 
and Fock shift terms (not shown here) also become 
weakly inhomogeneous around the vortices. 

7r-triplet Order Parameter: In the presence of coex- 



isting dSC and AFM orders, it had been argued [|32 34 
that an additional order, called 7r-triplet order, will de- 
velop self-consistently in the absence of any external mag- 
netic field. It was shown using mean-field gap-equations 
that even in the absence of any interaction that generates 
7r-triplet pairing directly, this order will be induced by co- 
existing DSC and AFM order. We define local 7r-triplet 
order as: 



H 



J 



{{ci+S[Ci]) ~ [cilCi+S])} e 



iQr 4 



(14) 



This defines a triplet pairing amplitude in the S z = 
channel. In the uniform system, Fourier transformation 
of Ilf results in (c_k^Ck+Qj), so that the center of mass 
momentum for the pair is Q = and hence the 

name. (The prefactor J/4 is arbitrary in the absence of 
any microscopic interaction that generates this order, we 
use this prefactor to compare the strength of this order 
relative to A*.) One physical way of understanding the 
induction of the 7r-triplet order is the folowing. The dSC 
quasiparticles carry momenta (k "f, -k j), whereas the 
AFM quasiparticles carry momenta (k + Q f, k f). As 
a result, in the presence of both dSC and AFM order, 
there is a possibility of pairing between quasiparticles 
with momenta (k + Q j, -k |) which leads to the tt- 
triplet pairing. 

Our t— J Hamiltonian does not contain any interaction 
involving at the mean field level (Eq. @), but the BdG 
self-consistency induces Ilf around the vortex, as shown 
in Fig. (|5|). On the other hand, the extended Hubbard 
model supports 7r-triplet pairing, and such an interaction 
should be present in the microscopic mean-field Hamil- 
tonian equivalent to Eq. (||). We have verified that the 
qualitative results from the extended Hubbard model are 
similar to the results presented here. 
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FIG. 5. Upper panel: 7r-triplet pairing amplitude IT; in- 
duced on a magnetic unit cell of size N — 26 x 26 system, with 
J = 1.15t and (n) — 0.875. 11; is zero at the vortex core, rises 
around the vortex, and falls gradually at larger distances. 
Lower panel: The grey-scale plot for 7r-triplet showing the 
length scale of decay which is similar to the AFM order (in 
Fig. |^). The dark (light) regions indicate larger (smaller) 
values for IT;. 11; also goes to zero on the lines of zero super- 
current. 

This Ilf order parameter was first introduced by Zhang 
p(| in a consistent picture for AFM and dSC order for 
cuprates using SO (5) theory. It was shown using general 
symmetry arguments that a theory which allows Af , Ilf* 
and rrii orders will satisfy the relation: 



Afnf = -mi(l - (m)), 



(15) 



where the variables with 'tilde' are the same as in Eq. (||) 
and in Eq. (|l4|) , but without the prefactors involving the 
coupling strength. This relation can be understood using 
a simpler picture of a so-called "spin-flop" model of an 
easy-axis AFM in a parallel magnetic field. Zhang |2(| 
discussed the analogy between his SO (5) model and the 
spin-flop model, where for small fields the equilibrium 



state is one with an AFM order parameter, N z , along 
the easy axis. As the field is increased above the spin- 
flop transition, the Neel vector flops into the plane with 
a value N± and there is also a uniform magnetization, 
M z , along the field. The case in which dSC and AFM 
order coexist corresponds to one in which the Neel vector 
N = {Nj_, N z ) is tipped out of the plane and the uniform 
magnetization M = (M±,M Z ) is tipped away from the 
easy axis. M± is analogous to the II triplet order pa- 
rameter. Since N is the difference of the two sublattice 
magnetizations, while M is their sum, N ■ M = 0. Us- 
ing the correspondence between the two models, m «-> 
N z , A «-> N±, (1 - (rii)) = n H <-► M z , U <-> M±, this 
implies the relation of Eq. ([Toj). 

As seen from Fig. (|]), (n VO rtcx) ~ 0.975 at the vor- 
tex center and hence does not achieve a value of true 
half-filling. However, calculations with increasing system 
sizes indicate that with larger system sizes the value ap- 
proaches closer to (n,) = 1. Defining the local hole den- 
sity as ((jivortcx) — ( n i)), instead of (1 — (n*)), we found 
that Eq. (|l5|) is satisfied within 5% accuracy. Although 
it is not obvious that our model has the same symme- 
try as the spin-flop model, we nevertheless find that this 
non-trivial relationship (Eq.^Hj) holds rather well site by 
sitel 



B. Local Density of States (LDOS) at the vortex 
core 

Having understood the self-consistent spatial struc- 
tures of the relevant order parameters in an antiferro- 
magnetic d-wave vortex, we turn our attention to the lo- 
cal density of states around the core. Study of the LDOS 
provides information about the low energy quasiparticle 
structure which is responsible for most of the novel phe- 
nomena associated with the cuprates. Also the LDOS is 
proportional to the tunneling conductance that is mea- 
sured directly in STM experiments |7],[l9| and hence can 
provide a direct test for the validity of theoretical pre- 
dictions. As discussed above, a pure d-wave BCS mean 
field calculation j| produces a large, unobserved ZBCP 
in the LDOS at the core of a d-wave vortex lattice. The 
ZBCP can be understood as the effect of Andreev re- 
flection of d-wave quasiparticles within the vortex core. 
Andreev reflection is caused by the rotation of the inter- 
nal state of an excitation from particle-type to hole-type 
(or vice versa) J33| . For a d-wave order parameter, such 
reflections scramble the phase information of the order 
parameter, and, as a result, the coherence peaks of the 
dSC collapse onto a ZBCP in the LDOS at the vortex 
core. Though the destruction of coherence peaks will oc- 
cur for a wide range of angle for a reflecting surface, the 
effect is most drastic for the angle of ir/4. The obvious 
question is: Does the additional AFM order suppress the 
unphysical ZBCP at the vortex core? To address this is- 
sue, we present the LDOS at the vortex core as obtained 



7 



from our calculations. 

The single particle local density of states at T — is 
calculated from the expression: 



1 



N ( i ,") = Jf £ {l< Q (k)| 2 <5( W -K(k)) 

n,k,a 

+ |< Q (k)| 2 <5(c + K(k))} 



(16) 



where the individual delta functions are broadened with 
a width comparable to the average energy level spacing. 
The results are shown in Fig.(^|). 



o.oi 



0.008 



0.006 



0.004 - 



0.002 




FIG. 6. Normalized density of states N(ui) at the vortex 
core for J = 1.15t and (n) = 0.875. The competition of 
dSC and AFM order splits the zero bias peak which would 
be present in the absence of AFM order. The signature of a 
spin gap appears at large ui ~ 0.5. The inset shows the DOS 
for a uniform system in the absence of vortices (y-axis levels 
are not given for clarity). The parameters for the inset are 
A max = 0.36, /i = 0.35 and m = 0.1, consistent with BdG 
results. 

There are two important features in Fig.(^|) that re- 
quire further explanations, (i) The zero bias peak splits 
in two, and the coherence peaks are suppressed, (ii) An 
extra feature appears at large positive to which looks like 
a second gap. Also, N(u>) appears to be quite asymmet- 
ric. 

The asymmetry and the second gap scale (at lo ~ O.bt) 
in the LDOS can be understood by analyzing the problem 
of a uniform d-wave SC coexisting with AFM order in the 
absence of vortex. This is described by the Hamiltonian 
H = EL *k M k*k where, 



I £ k - n A k -to n k 

A k -£k + /i -n k -m 

-to -n k -f k - fi -A k 

V n k to -A k & + /■«/ 



(17) 



*k = (ckt, cL ki , c k+QT , c^ k _ QT ^ and the prime on 
the k-sum indicates a sum over only half of the Brillouin 



zone, which has been reduced by AFM order. Here £ k = 
—2t(cosk x + cosky), A k = At^cosfca, — cosk y ) and Il k = 
Ilc^cosfca; — cosky). The result of the spin-density wave 
(SDW) problem with A = 0, n = and fi = is well 
known; under these conditions diagonalization of Eq. jl^ ) 
results in a N(u>) that has a spin gap of size m at the 
Fermi energy. Inclusion of particle hole asymmetry (/j, ^ 
0) shifts the overall density of states by an energy, /i. 
Superconductivity can now be introduced with a non- 
zero Ao, which will open a superconducting gap at lo = 
with d-wave symmetry. The resulting DOS for this 
uniform system, which is plotted as an inset of Fig. (|^), 
shows two gap scales, one at to — and the other at lo = 
/i| +m. With to = 0, a non-zero IIo in Eq. ( p7| ) splits the 
dSC coherence peaks at A max ±IIo; the splitting becomes 
asymmetric when fi ^ 0. However, the 7r-triplet order 
does not change N(us) qualitatively from that in the inset 
of Fig. (||) (at least for the choice of our parameters). (See 
also Fig. (10) of Ref. The presence of the vortex 

modifies the LDOS structure near to = (described in the 
next paragraph), but the SDW gap at larger to survives 
the vortex excitations and shows up in Fig. (||). The 
above argumnet for the origin of the 'spin-gap scale' was 
verified by tuning model parameters that change /i and 
to. It was found that the 'spin-gap' changes accordingly 
in our BdG solutions. 

Next we turn our attention to the splitting of the zero 
bias peak in N(io). We have already argued (See Sec. II) 
that the ZBCP in the LDOS in Ref. {§ is due to the An- 
dreev reflection of d-wave quasiparticles from the vortex 
core which acts like a boundary. The effects of an addi- 
tional subdominant ordering on the ZBCP at the surface 
have been investigated recently 36 - 39 . The conclusion 
is that, if the subdominant order breaks time-reversal 
symmetry, then the ZBCP will split in two. For our 
case, the subdominant order is AFM order, which breaks 
time-reversal symmetry, so that the presence of the AFM 
order within vortex core should split the ZBCP, as found 
in Fig. (^). Within this formalism, the energy differ- 
ence between the split peaks would be proportional to 
the value of rrii at the vortex core, which is consistent 
with our findings. 

Having understood the features of the LDOS at the 
vortex core as found in our calculations, we must em- 
phasize that, a mean-field calculation on a model with 
static antiferromagnetism like ours, does not reproduce 
all the experimental findings. Although the competition 
between dSC and AFM order produces a gap near the 
Fermi energy in the DOS - a result consistent with ex- 
perimental findings - the spin-gap at larger lo has not 
been observed in experiments. Also, the strong split- 
peak around lo = 0, which is due to static AFM does not 
resemble the weak subgap hump found in tunneling con- 
ductance data. The strong split peaks, however, die off 
rapidly away from the vortex center, as seen from the top 
panel of Fig. (^|), which shows the evolution of the LDOS 
along the diagonal direction from the vortex. We also 
find that the coherence peaks build up near lo ~ A max 
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away from the core, as expected. Similar behaviour is 
found along the x and y directions as well. 



N(w) - 



- tj 

"r TOrt + (1.1 ),_,.; \ >\ n t 

r yort + (3,3) . /V\J 


- 









0.004 



0.003 



DOS averaged within 
vortex core region _ 




FIG. 7. Upper panel: Local Density of states away from 
the vortex core along the diagonal direction. The strong 
split-peaks fall rapidly in magnitude away from the vortex 
center. 

Lower panel: Normalized local density of states averaged over 
all sites within a radius £ from the vortex center. 

Surprisingly the low energy peaks persist to larger dis- 
tances, even outside the vortex core. This is interesting 
from the experimental point of view, because such low 
energy quasiparticle humps have been observed beyond 
the vortex core fisf] . This led us to calculate the local 
density of states averaged over all sites within a radius 
£ from the vortex core, which is presented in the lower 
panel of Fig. (^) . Allowing for some uncertainty in locat- 
ing the vortex center in STM measurements, we believe 
the qualitative agreement of our results in Fig. (]?]) with 
experimental findings is significant. 



IV. CONCLUSIONS 

In conclusion, we have studied the spatial structures 
of different order parameters and their effect on the local 
density of states at the vortex core within a model that 
allows for the competition between d-wave superconduc- 



tivity and static antiferromagnetism. Our main results 
are: (1) Antiferromagnetic order develops near the vortex 
center, where the d-wave order parameter is suppressed, 
and exhibits interesting structure within the magnetic 
unit cell. (2) The length scale for decay of AFM order 
is larger than the superconducting coherence length. (3) 
Coexistence of dSC + AFM order induces TT-triplet order 
in the vortex core region. (4) Competition between dSC 
and AFM splits the zero bias peak in the tunneling con- 
ductance. Experimental data strongly suggest that, at 
the vortex core, the low energy states near gap nodes are 
moved to higher energies. We believe that such physics 
can not be fully captured within mean-field theory. Dy- 
namic AFM fluctuations and strong correlation effects 
should be included to make a more detailed comparison 
to the experimental data. 

Appendix 

Self-consistent Parameters at the Boundary 



The boundary terms corresponding to Eq. (10) will be 
modified by making the following replacement on the left 
hand side of Eq. ( |10| ) , 



L i±y^N y y 



a(k)e 



(18) 



and the definition of A^ n (i) will be modified in the 
boundary as following: 



A ±x 
l qn 



(0 



Je 



±iai v N x 



AN 



[{ 



^ik^N^^n 



«u(q)«?±*FJv.*,t(q)} (i - /qni) 

+ {e Ti ^^< T (q)<± £T7 v x , a (q) 
^^^ T ( q )<t(q)}/knT] 

-{ e T^ <± . T ^, a(qK * (q) 

'<i(q)«,W 9 »,t(q)} (i - /q«i) 



e 
J 

w 

±iq y N v n 



{e^^»< T (q)C iTWtW (q) 
e ±^iV Vn » T(q)u n* (k) | /qnT j 



(19) 



Self-consistency conditions for Hartree and Fock shifts 
(equivalent for Eq. (|l2])) are not given here; they are 
easily derivable in our formalism. 
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